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An improved version of the PNJL model is used to calculate various thermodynamical 
quantities, viz., quark number susceptibility, isospin susceptibility, specific heat, speed of 
sound and conformal measure. Comparison with Lattice data is found to be encouraging. 



I. INTRODUCTION 

Under the extreme conditions of temperature and/or density a phase transition is expected in 
the QCD phase diagram. In fact, there are two phase transitions at these extreme conditions, 
namely the confinement to deconfinemcnt phase transition and the chiral phase transition. These 
two phase transitions are defined in two extreme quark mass limits. For infinite quark mass 
limit confinement to deconfinement phase transition is well defined, and Polyakov loop acts as the 
order parameter of the transition corresponding to the Z{3) (for three colours) global symmetry 
breaking. In the limit of zero quark mass we have chiral symmetry restoring transition. This phase 
transition is characterized by the chiral condensate, which acts as the order parameter. But, in 
the real world where quark masses have a finite value, the nature of these phase transitions is not 
known. Whether these two transitions will occur simultaneously or one will precedes the other 
remain an open question. Order of the phase is another aspect of the phase transition which is 
also being pursued actively. 

Here we are concerned with the thermodynamic aspect of the phase transition in terms of 
some relevant thermodynamic variables such as various quark number and isospin susceptibilities, 
specific heat, speed of sound etc., to see how the behaviour of these variables change as we go 
from one phase to the another. Being in the nonperturbative domain the Lattice QCD and the 
effective models are the only theoretical framework which can be employed to study the phase 
transition. The effective model we employ here is known as the PNJL (Polyakov loop + Nambu- 
Jona-Lasinio) model. Within the framework of PNJL model one can study both the confinement 
to deconfinement transition and the chiral phase transition together. For details of this model and 
of its various improvements see Refs. 1, 2, 3, 4, 5, 6, 7, 8, 

The plan of the paper is as follows. In the next section we briefly outline the theoretical frame- 
work. In section 3, our results are discussed and finally we conclude in section 4. 
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II. THEORETICAL FRAMEWORK 



The coupling between the chiral and the deconfinement order parameter within the PNJL model 
allows us to study the thermodynamics of both the transitions inside a single theoretical framework. 
The thermodynamics of the gauge sector is controlled by the Polyakov loop whereas the NJL 
model takes care of the thermodynamics of the quark sector. To obtain various thermodynamic 
quantities we start with the thermodynamic potential per unit volume computed within the PNJL 
model : 
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where au and ad are, respectively, two light flavours condensates and the respective chemical 
potentials are /i^ and Note that ~ + Md)/2 and — — /^(i)/2. The quasi-particle 

energies are Eu,d = \Jp ^ + ^, where mu,d — wo — 4:Giau,d — ^G2(Jd,u are the constituent 

quark masses and mo is the current quark mass (we assume flavour degeneracy). Gi and G2 
are the effective coupling strengths of a local, chiral symmetric four-point interaction. We take 
Gi = G2 = G/4, where G is the coupling used in Ref. '3"|. A is the 3-momentum cutoff in the NJL 
model. The form of the effective gauge potential is in the form given by, 

iY($,$,T) W (*,$,T) , , 

y4 - r4 -^ln[J($,1')] , (2) 

where W is taken from [3i, iSi] . It arises because of the transformation from matrix valued field L to 
complex valued field 4>. The role of this second term (known as Vandermonde determinant (VdM) 
term) is to regulate the behaviour of $ so that its value remains within the group theoretic limit. 
This term is dropped while calculating the pressure. For details we refer to [9|. 

Once the thermodynamic potential is fixed, we first obtain the mean fields by minimizing 17 
w.r.t the field variables <&, au and ad- We then calculate all the relevant thermodynamic quan- 
tities by appropriate differentiation of fi. Here we are concerned with the various susceptibilities 
namely the quark number susceptibility (QNS) and isospin number susceptibility (INS). It is also 
to be noted that by QGP susceptibility we intend to mean these susceptibilities. Susceptibilities 
are important as they are related to fluctuations which can be measured experimentally. 

We computed all these quantities from the Taylor expansion of the pressure, as it is related to 
the thermodynamical potential via the relation, 

P{T,ii^) = -n[T,iio) ■ (3) 

For susceptibilities we expand the scaled pressure (P/T*) in a Taylor series for the quark number 
and isospin number chemical potentials, fj,o and fj,i respectively. The coefficients we are interested 
in are given by. 



n 



I 9" {P{T,iio)/T^) 
I 



„nO 



dm 



(4) 



Mo=0 



, ^ 1 9" (p(r,/.o,M/)/r") 



- „(n-2)2 



; n > 1. 



(5) 



/iO=0,/J7=0 



The first equation represents the QNS and its higher order derivatives whereas the second one 
corresponds to the INS and its higher order derivatives. Specific heat and the speed of sound are 
calculated from the expansion of the pressure with respect to the temperature. Explicit expressions 
of these quantities are given by: 
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For details of computing these coefficients please refer to [1, ■ 



III. RESULT 



(7) 



1 








1 0.35 










0.3 


0.8 






H 3 B 


0.25 


0.6 








0.2 


0.4 






Co (Lattice) 


0.15 


0.2 








Cj (Lattice) 
cf(PNJL) — 
cf (PNJL) 


0.1 
0.05 





0.6 0.8 


1 1.2 1.4 
T/Tc 


1.6 1.8 2 



C4 (Lattice) i-e-i 
C4 (Lattice) b 
cf (PNJL) 




(a) 



(b) 



FIG. 1: The Taylor expansion coefficients of the pressure in quark number and isospin chemical potentials 
as functions of T/Tc. Symbols are LQCD data Arrows on the right indicate the corresponding ideal 
gas values. 



We present the QNS, INS and their higher order derivatives with respect to /iq in Fig. [T] We 
have plotted the LQCD data from Ref. 12] for quantitative comparison. We find, QNS C2 (Fig. 



1(a) I matches well with the Lattice data upto \.2Tc but beyond which it goes to SB limit whereas 



the Lattice data saturates at 80% of the ideal gas value. Similar trend is observed for INS C2. It 
remains close to LQCD data upto 1.2Tc and then it saturates at SB limit. However both QNS C2 
and INS C2 agree well with each other which is also observed in Lattice calculation. It is to be 
noted here the QNS C2 matches well with the LQCD data as shown in [l3|, if we do not include 
the VDM term. 

Variation of next higher order coefficients C4 of both QNS and INS with temperature are shown 



in Fig. 1(b) We find a good agreement for both the coefficients with the Lattice data for whole 



range of temperatures. If we do not include the VDM term (see (lOj) then the QNS C4 matches the 
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FIG. 2: Temperature dependence of (a): energy density e and specific heat Cv- (b): squared speed of 
sound v1 and conformal measure A/e. The arrows in the right show the corresponding SB hmit. 

Lattice data only upto T ^ l.lTc- The effect of the VDM term is to bring the C2 and C4 coefficients 
to their respective SB hmit. 

In Fig. [5] we have presented our resuh for specific heat, speed of sound and conformal measure. 
Compared to our previous observation in Ref. the specific heat Cv differs only marginally at 
higher temperature. We find a small peak near the transition temperature which in the case of 
critical end point (CEP) will diverge. At higher temperature it saturates below the ideal gas limit 
and approaches the conformal limit of Ae/T'^ from above. The effect of the VDM term cancels out 
for both the speed of sound and the conformal measure as they are the ratios of the energy density 
and pressure. 

IV. SUMMARY 

Various thermodynamical variables of importance are calculated using an improved version of 
the PNJL model. All the variables estimated are found to be in reasonable agreement with the 
LQCD data, though the QNS C2 and the INS on the Lattice are smaller by about 20 %. We found 
all thermodynamical variables approach to their respective ideal gas limit at higher temperature. 
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